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Abstract. This paper presents infrared data obtained from observations carried out at the ESO 3.5m New Technology Telescope 
(NTT) of the Hubble Deep Field South (HDF-S) and the Chandra Deep Field South (CDF-S). These data were taken as part of 
the ESO Imaging Survey (EIS) program, a public survey conducted by ESO to promote follow-up observations with the VLT. In 
the HDF-S field the infrared observations cover an area of ~ 53 square arcmin, encompassing the HST WFPC2 and STIS fields, 
in the JHK S passbands. The seeing measured in the final stacked images ranges from 0"79 to L'22 and the median limiting 
magnitudes (AB system, 2" aperture, 5 G detection limit) are Jab ~ 23.0, Hab ~ 22.8 and Kab ~ 23.0 mag. Less complete 
data are also available in JK S for the adjacent HST NICMOS field. For CDF-S, the infrared observations cover a total area of 
~ 100 square arcmin, reaching median limiting magnitudes (as defined above) of Jab ~ 23.6 and Kab ~ 22.7 mag. For one 
CDF-S field H band data are also available. This paper describes the observations and presents the results of new reductions 
carried out entirely through the un-supervised, high-throughput EIS Data Reduction System and its associated EIS/MVM C++- 
based image processing library developed, over the past 5 years, by the EIS project and now publicly available. The paper also 
presents source catalogs extracted from the final co-added images which are used to evaluate the scientific quality of the survey 
products, and hence the performance of the software. This is done comparing the results obtained in the present work with those 
obtained by other authors from independent data and/or reductions carried out with different software packages and techniques. 
The final science-grade catalogs together with the astrometrically and photometrically calibrated co-added images are available 
at CDS**. 
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1. Introduction 

One of the goals of the ESO Imaging Survey (EIS, 
Renzini & da Cost al ll997l) was to carry out moderately deep, 
multi-passband, optical/infrared observations over relatively 
large areas to produce faint galaxy samples. The primary ob- 
jective was to produce samples with color information to esti- 
mate photometric redshifts, extract statistical samples of galax- 
ies likely to be in the poorly sampled 1 < z < 2 redshift interval 



Send offprint requests to: L.F. Olsen, LisbethFogh.01sen@obs-nice.fr 
* Based on observations carried out at the European Southern 
Observatory, La Silla, Chile under program Nos. 61.E-9005, 162. 0- 
0917, 163.O-0740 

** The science grade images and catalogs are available at CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ AJ 



or Lyman-break candidates at z > 2.5 and, perhaps, examples 
of high redshift quasars and galaxies all of which are interest- 
ing targets for follow-up spectroscopic observations with the 
VLT. 

The first such attempt was the EIS-DEEP survey carried out 
at the NTT in the period 1998-2000. The infrared part of the 
survey was carried out in the period 1998-1999 and is reported 
in the present paper. While by recent standards this survey is 
relatively modest in size, the combination of depth and area 
made it unique at the time it was conceived. The survey con- 
sisted of optical and infrared observat i ons o f the Hubble Deep 
Field South (HDF-S, IWilliams et aD. l2000[) and the Chandra 
Deep Field South (CDF-S. iGiacconi et alll200lh . The choice 
of the regions was an attempt by the EIS Working Group to re- 
concile the great interest generated by the HST observations of 
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a region selected in the southern hemisphere (HDF-S), with the 
desire to conduct observations in a less crowded region and, in 
particular, devoid of bright stars. Therefore, the second region 
was chosen to overlap the area, selected for its low HI column 
density, for which deep Chandra observations were already be- 
ing planned. This region is particularly well-suited for ground- 
based imaging due to its lack of bright stars. Not surprisingly, it 
eventually became a widely explored window to the high red- 
shift Universe with deep multi-wavelength observations being 
conducted by all major space obse rvatories and several g round- 
based follow-up programs (e.g. iGiavalisco et alll2004l) 

Preliminary r esults of this su r vey w ere original ly reported 
in pre prints by Ida Costa et alJ Jl998l) and iRengelink et alJ 
(1998) who presented fully calibrated optical and infrared im- 
ages, single-passband catalogs extracted from these images, 
color catalogs based on a reference % 2 image and lists of color- 
selected high-redshift galaxy candidates utilizing the Ly-break 
technique. All of these data products were made publicly avail- 
able world-wide. The optical images were reduced using IRAF 
and the infrar ed images using t he program Jitter of the Eclipse 
library (Devillard et ali Il999h developed by ESO to reduce 
ISAAC and SOFI data. 

Despite satisfactory attempts to assess the quality of the de- 
rived products by, for instance, comparing counts of galaxies 
with those of other authors and of star counts and the stellar lo- 
cus in color-color diagrams with model predictions, subsequent 
comparisons with deeper data and with independent analysis 
using other reduction techniques showed that the reductions 
carried out by the EIS team using the available Jitter package 
led to a systematic loss of flux at faint magnitudes. The color 
catalogs released also showed some unexpected results which, 
after some time, were traced to relative offsets between the as- 
trometric calibration of the optical and infrared images. Due to 
these problems the abovementioned preprints were withdrawn 
and the data release described in the present work supersedes 
those made earlier by the EIS team of this dataset. The iden- 
tified problems, together with the need to develop a common, 
high-throughput system to handle optical/infrared, single and 
multi-chip instruments, led to the development, as part of the 
EIS proje ct, of a n image processing engine (EIS/MVM library, 
Vandame, 2004). This library is a central part of the integrated , 
end-to-end EIS Data Reduction System Jda Costa et all l2004) 
built to carry out un-supervised image reductions, stacking and 
catalog preparation. Preliminary results obtained with the new 
image processing library, in particular addressi ng the flux los s 
problem were originally reported by Vandam e et alJ l|2001). 
However, even though a significant improvement was achieved, 
additional testing indicated that an error was introduced in the 
code. To avoid additional intermediate releases, the final re- 
lease of these data was postponed until the completion of the 
EIS Data Reduction System infrastructure which includes the 
photometric pipeline and provides the infrastructure required 
for re-processing and re-calibrating data in an efficient way, a 
key element for testing new software. 

In this paper these new reductions of the data accumulated 
by EIS-DEEP, including new ones not previously released, are 
presented. In Sect. |2] the observations, including those carried 
out in 1999 and not previously reported, are reviewed. The goal 



Table 1. SOFI Pointings (J2000.0). 



Field 


a 


8 


HDF-S- 1 


22:33:29.1 


-60:33:50 


HDF-S-2 


22:32:42.4 


-60:33:50 


HDF-S-3 


22:33:00.0 


-60:37:59 


CDF-S- 1 


03:32:16.7 


-27:46:00 


CDF-S-2 


03:32:38.0 


-27:46:00 


CDF-S-3 


03:32:16.7 


-27:50:35 


CDF-S-4 


03:32:38.0 


-27:50:35 



Table 2. Number of observed fields in each region. 



Region 


J 


H 


K s 


CDF-S 


4 


1 


4 


HDF-S 


3 


2 


3 



of this paper is two-fold: first, provide the most recent reduc- 
tions and substitute previous releases; second, use these data as 
a test case to justify the use of the developed software to treat 
the much larger dataset accumulated by the infrared part of the 
Deep Public Survey (DPS) carried out by the EIS team over 81 
nights in the period 2000-2004, which generated about 15000 
science frames. The results for thi s survey will be pre sented in a 
forthcoming paper of this series dOlsen etall E006). In Sect. [3] 
the procedures adopted to remove the instrumental signature 
and to astrometrically and photometrically calibrate the nightly 
data are discussed, while in Sect.|4]the attributes of the final co- 
added images and their associated catalogs are presented. The 
scientific quality of these advanced survey products is assessed 
by comparing in Sect. |5] the results derived from the present 
dataset with those obtained from a variety of other datasets re- 
duced using diverse techniques. The main results of the present 
paper are summarized in Sect. [6] 

2. Observations 

As described above, the EIS-DEEP observations originally 
were planned to cover two regions: one region consisting of 
three adjacent fields, of about 25 square arcmin each, covering 
the WFPC2, STIS and NIC3 fields; and another region consti- 
tuting a 2 x 2 mosaic, comprising a total area of ~ 100 square 
arcmin covering part of the now well-known CDF-S region. 
The center of the mosaic was chosen to coincide with the nomi- 
nal position selected for the X-ray observations, where the sen- 
sitivity and image quality were expected to be the best. The 
pointings adopted in the construction of the mosaics are listed 
in TableG] 

The observations were obtained using the SOFI cam- 
era JMoorwood et all Il998l) mounted on the NTT. SOFI is 
equipped with a Rockwell 1024 2 detector that, when used to- 
gether with its large field objective, provides images with a 
pixel scale of 0.29 arcsec, and a field of view of about 4.9 x 
4.9 square arcmin. The total area covered in J and Ks is about 
~ 175 square arcmin split into two mosaics, one of ~ 75 square 
arcmin covering HDF-S and the other of ~ 100 square arc- 
min covering CDF-S. The observations were carried out in 
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groups of exposures (OBs). The excution of an OB consists 
of a series of short exposures with small position offsets from 
the target position. The observations were conducted using the 
Auto Jitter observing template, which implements the off- 
setting of adjacent exposures within a box of a user specified 
size s, chosen to be s = 45 arcsec, approximately 15% of the 
field of view. The offsets are constrained so that all distances 
between pointings, in a series of 15 consecutive exposures are 
larger than 9 arcsec. Each observation comprised sixty one- 
minute exposures, which themselves were the average of six 
ten-second sub-exposures. Table [2] summarizes the available 
observations for each of the mosaics. The numbers listed in 
the table represent the number of observed fields regardless of 
their exposure times. 

A total of 3051 science frames were obtained. Out of 
those, 1316 are in the CDF-S region and represent 22 hours 
of on-target integration, obtained during 7 nights in the period 
August 3, 1998 and November 11, 1998. The remaining 1735 
are located in the HDF-S region and were accumulated in 29 
hours of on-target integration, in observations carried out dur- 
ing 14 nights over the period August 3, 1998 to July 7, 1999. 
Observations of standard stars are available for all nights, gen- 
erating 613 images taken in 7805 sec. 

3. Data reduction 

The new reduction of the data, described in the present pa- 
per, was carried out using the EIS Data Reduction System 
(da Costa et al ., 2004), which is an end-to-end, integrated sys- 
tem providing the infrastructure for the un-supervised reduc- 
tion and calibration of large volumes of imaging data. An in- 
tegral part of the reduction system is a C++-based image pro- 
cessing library (EIS/MVM library, Vandame, 2004) developed 
as part of the EIS effort. It is both a library of basic image pro- 
cessing routines, and an integrated, high-throughput pipeline 
(reaching a sustainable rate of about 0.5 Mpix/sec) for auto- 
mated reduction of single or multi-chip optical/infrared images. 

The main steps of the EIS Data Reduction System are: the 
creation of reduced images on a night by night basis; the pho- 
tometric calibration of the reduced images; stacking of the re- 
duced images yielding the final stacked images forming the 
basis of the catalog extraction. The EIS/MVM pipeline is re- 
sponsible for the first step of producing fully reduced im- 
ages and weight-maps carrying out bias subtraction, flatfield- 
ing, de-fringing, correction of eletronic cross-talk effects, iden- 
tification and automatic masking of satellite tracks, and sky- 
background subtraction as well as first-order pixel-based im- 
age stacking (allowing for translation, rotation and stretching of 
the image). The stacking process combines the reduced images 
into final image stacks properly flux scaling the images to con- 
serve the photometric calibration. Finally, the object catalogs 
are extracted using SExtractor. The reducti ons of optical wide- 
field imaging data are present ed by Dietrich et al 1 il2005[ here- 
after Paper I) and Mi gnano et alJ ( 120061) . where details about 
common aspects of the reduction can be found. Paper I also 
includes a detailed description of the released images and cat- 
alogs and their format. Here only the more important steps for 
handling infrared data are reviewed. 



After the basic image reduction, the de-fringing, sky- 
subtraction and elimination of image artifacts are carried out in 
a two-step procedure: first, a preliminary image is obtained by 
stacking a set of sky-subtracted frames, with the background 
estimated using a running median; second, objects are iden- 
tified in this stacked image and masks are created and en- 
larged by a user-specified factor; finally, the de-fringing/sky- 
subtraction step is repeated with the object masks mapped to 
the jittered position of each individual image. In the final co- 
addition of the images the masked pixels are neglected and 
deviant pixels are clipped. This procedure is particularly im- 
portant for infrared observations since, as explained in Sect. [2] 
the observations are carried out as sequences of short expo- 
sures, in order to allow a suitable estimate of the fast varying 
background at these wavelengths to be made. The short expo- 
sure time of the individual images, prevent the detection of a 
large number of objects with non-negligible flux, which can 
make a significant contribution to the background when several 
images are co-added. Therefore, if these objects are not prop- 
erly masked when estimating the background, the background 
is over-estimated, biasing the flux of even relatively bright 
objects. This was the effect observed in the original release 
of the EIS-DEEP data using the Eclipse red uction package 
(Ida Costa et all Il998l iRengelinket ail Il998b . With the two- 
step procedure now adopted, the flux estimates have greatly 
improved and, as shown in Sect. [5] the flux estimates are in 
excellent agreement with those obtained using other reduction 
packages commonly used to reduce infrared data. 

The library also includes advanced algorithms used for the 
astrometric calibration and warping of reduced images onto a 
user-defined reference grid. The astrometric calibration of the 
infrared frames was performed us ing the GSC2.2 (HDF-S) and 
USNO-B (CDF-S, Monet et all 120031) catalogs as refe rence 
following the method developed by Diam diTet all {l993). The 
GSC2 catalog was used because of the lack of reference stars in 
the USNO-B within the surveyed are of the HDF-S region. The 
internal astrometric accuracy is about half a pixel (^0.15 arc- 
sec) and is limited by the internal accuracy of the reference cat- 
alog. When dealing with multi-wavelength surveys it is prefer- 
able to register the infrared data using a catalog from an optical 
counterpart image. This will not only provide a much larger 
number of reference objects but also avoid the effects associ- 
ated with small deviations of the astrometric solutions obtained 
for the optical and infrared data that were identified in the orig- 
inal release of the EIS-DEEP data. 

The photometric calibration for each night is based on ob - 
servations of standard stars taken from IPersson et all (^998). 
Typically two to seven stars were observed over a range of air- 
mass. For all nights, independent photometric solutions were 
attempted using the photometric pipeline of the EIS Data 
Reduction System. Depending on the airmass and color cov- 
erage linear fits with one, two or three free parameters are at- 
tempted. If no standard star observations are available or no 
suitable fits are possible a default solution is adopted. This was 
the case for only two nights: August 4, 1998 for both / and 
K s bands and November 11, 1998 for the /^-band with insuffi- 
cient observations of standard stars during the night to derive a 
reliable solution. 
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Table 3. Type of photometric solutions available for each com- 
bination of mosaic and band. 



Region Passband Default 1-par 2-par 3-par total 



Table 6. Quality assessment of reduced images. 



CDF-S 
CDF-S 
CDF-S 
HDF-S 
HDF-S 
HDF-S 



./ 
H 

K S 
./ 
H 



Table 4. Average photometric solutions. There was no case 
where the color term could be determined. 



Passband 


Zp 


k 


J 


23.21 ±0.05 


0.11 


H 


23.00 ±0.05 


0.02 


K s 


22.45 ±0.06 


0.13 



Table 5. Comparison between EIS and Telescope Team solu- 
tions (EIS -Telescope). There was no case where the color term 
could be determined. 



Passband 


Zp 


k 


J 


-0.07 


0.01 


H 


-0.06 


-0.02 


K s 


0.0 


0.03 



Table[5]summarizes the type of photometric solutions avail- 
able for each mosaic. The table gives in Col. 1 the region name; 
in Col. 2 the passband; in Col. 3 the number of default solu- 
tions; in Cols. 4-6 the number of nights with photometric so- 
lutions derived from 1- to 3- parameter fits; and in Col. 7 the 
total number of nights with science and standard star observa- 
tions for a given combination of mosaic and passband. It can 
be seen that in many cases only 1 -parameter fits are available, 
in particular in the HDF-S mosaic, and that no solution was ob- 
tained with an independent estimate of the color term. This is 
the result of insufficient airmass and/or color coverage of the 
standard star observations. Note that for mosaics, in principle, 
one good solution suffices to determine a reliable calibration, 
using overlapping regions to propagate a given solution. 

Table 0] summarizes the average photometric solutions de- 
termined for the entire period of observations, covering more 
than one year. The table lists in Col. 1 the passband; in Col. 2 
the mean zeropoint and standard deviation; in Col. 3 the mean 
extinction. The number of solutions is in general small rang- 
ing from 1 to 6 solutions in the different bands and only for 
1-3 of these could the extinction be determined, as can also 
be seen from Table [3] Therefore, the standard deviation for the 
extinction values could not be determined reliably. The errors 
estimated for the extinction values are typically of the order 
of < 0.05. The listed standard deviations are a good indication 
of the uncertainty of the estimated zeropoints. The solutions 
have been compared with the values provided by the Telescope 



Region 


Passband 


A 


B 


C 


D 


CDF-S 


J 


13 











CDF-S 


H 


1 











CDF-S 


K s 


10 


3 








HDF-S 


.1 


11 











HDF-S 


H 


7 


2 








HDF-S 


K s 


3 


12 









Team 1 in Table [5] listing in Col. 1 the passband; in Col. 2 the 
difference in zeropoint of the photometric solutions computed 
as (EIS-Telescope Team); in Col. 3 the corresponding differ- 
ence for the determined extinction. It can be seen that the dif- 
ferences are all within the estimated uncertainty, which is note- 
worthy given that the solutions are not derived for the same pe- 
riod of time. Below, the derived magnitudes are compared with 
those of other authors, and thus careful consideration is given 
to the photometric quality of the data. Where the magnitudes 
are listed in the AB system the following conversions have been 
used: J AB = J ±0.904, H AB =H+l .374 and K AB = K s + 1.841 . 

Before creating stacks and associated catalogs, discussed 
in the next section, all the reduced images were inspected and 
graded. Table[6]summarizes the results of the visual inspection 
grouping by mosaic. The table lists in Col. 1 the region name; 
in Col. 2 the passband; and in Cols. 3-6 the grades given rang- 
ing from A (best) to D (worst). These grades reflect the overall 
cosmetic quality of the images including the background and 
the presence of other features that may affect the source extrac- 
tion. In general, a grade A image shows no significant cosmetic 
problems, while grade D indicates that this image is useless 
for any further use. Grades B and C are intermediate, where C 
is usually not useful in itself but can still add information in 
the subsequent stacking procedure. In addition to the grade a 
subjective comment, to describe any abnormality or feature, is 
normally associated with the image. As can be seen from the 
table all images were graded A or B, corresponding to a high 
cosmetic quality. 

The 3051 science frames were input to the EIS Data 
Reduction System and converted into 62 reduced images, each 
corresponding to an observing block, with 35 in the HDF-S and 
27 in the CDF-S regions. Table summarizes the information 
on the available reduced images. The table gives in Col. 1 the 
region name; in Col. 2 the passband; in Col. 3 the number of 
fields observed; in Col. 4 the total on-source time in ksec; in 
Col. 5 the number of nights in which the observations were 
carried out; in Col. 6 the number of raw exposures involved; in 
Col. 7 the number of reduced images available. 

Typically, the reduction data rate is of the order of 0.2 
Mpix/sec for each dual processor. Since the system was oper- 
ating with 8 independent machines the effective data rate was 
of the order of 1.6 Mpix/sec. The high-throughput of the sys- 
tem together with its administrative infrastructure provided by 
its associated database and system architecture was essential, 
allowing re-processing and re-calibration, in particular, taking 
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Table 7. Summary of reductions. 



Region Passband # Fields T (ksec) Nights Exposures # Images 



CDF-S 


J 


4 


40.1 


3 


668 


13 


CDF-S 


H 


1 


3.5 


1 


58 


1 


CDF-S 


K s 


4 


35.4 


4 


590 


13 


HDF-S 


.1 


3 


29.9 


5 


498 


11 


HDF-S 


H 


2 


23.2 


4 


387 


9 


HDF-S 


K s 


3 


51.0 


6 


850 


15 



into account the large number of indivdual programs carried 
out by the EIS team. This, together with the infrastructure avail- 
able for keeping different versions of products, enabled the 
comparison of different reductions in an easy way, a key ele- 
ment for long term programs. 



4. Final products 

4. 1. Images 

The 62 reduced images mentioned in the previous section were 
converted into 17 stacked (co-added) images following the pro- 
cedures described in Paper I. All stacked images were visually 
inspected and graded. Again the quality was found to be good 
with 14 images being graded A and the remaining 3 B. The re- 
sulting mosaics in each passband are shown in Figs.^and[2]for 
HDF-S and CDF-S, respectively. 

The main attributes of the stacks produced for each mosaic 
are summarized in Table [8] The table gives in Col. 1 the field 
name; in Col. 2 the passband; in Col. 3 the grade; in Cols. 4 and 
5 the point-spread function (PSF) FWHM in arcseconds and its 
anisotropy (PSF rms) as measured in the final stacked image; 
in Col. 6 the limiting magnitude, m/, m , estimated for the final 
image stack for a 2"aperture, 5a detection limit in the Vega sys- 
tem; in Col. 7 the completeness of the observations expressed 
as the fraction (in percentage) of observing time relative to that 
originally planned. It can be seen that the seeing is, in general, 
quite good, ranging from 0.64 to 1.22 arcsec with a median 
of 0.84 arcsec. The median limiting magnitudes in the Vega 
system are 22.72, 21.47 and 20.90 in the J, H, and A^-band, 
corresponding to 23.62, 22.84 and 22.74 in the AB-system. 

As a final check of the quality of the photometric cali- 
bration, the regions of overlap between adjacent frames were 
used to estimate the relative field-to-field magnitude variations, 
yielding an estimate of the overall uncertainty in the absolute 
calibration. In the HDF-S there are only two overlaps in J and 
K s , with the mean values being 0.05 and 0.07, and comparable 
scatter. In the //-band there is only one overlap with an offset 
of 0.05. For the CDF-S there are four overlaps for the J and 
K s data, yielding a mean offset of 0.03 and 0.02 mag and a 
scatter of 0.029 and 0.024 mag, respectively. From these re- 
sults one can estimate the accuracy of the absolute calibration 
of HDF-S is ~ 0.05 mag, while that of CDF-S is considerably 
better < 0.03 mag in both bands. 
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Fig. 1. Mosaic image of the HDF-S in the J, H and K s -band 
from top to bottom. 

4.2. Catalogs 

Catalogs were extracted from the stacked images using the EIS 
Data Reduction System adopting the same procedure as de- 
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Table 8. Main characteristics of the final stacks. 



Field 


Passband 


Grade 


PSF FWHM 


PSF rms 


™lim ( Ve g a ) 


Completeness 


CDF-S- 1 


J 


A 


0.84 


0.194 


22.72 


94% 


CDF-S- 1 


K s 


A 


0.97 


0.183 


21.08 


105% 


CDF-S-2 


J 


A 


0.92 


0.139 


22.86 


89% 


CDF-S-2 


K s 


A 


0.70 


0.189 


20.90 


90% 


CDF-S-3 


J 


A 


0.80 


0.194 


22.95 


99% 


CDF-S-3 


K s 


B 


1.12 


0.123 


20.75 


84% 


CDF-S-4 


.1 


A 


0.64 


0.111 


21.44 


89% 


CDF-S-4 


H 


A 


0.84 


0.113 


21.21 


32% 


CDF-S-4 


K s 


A 


0.87 


0.122 


20.64 


49% 


HDF-S- 1 


.1 


A 


1.22 


0.113 


22.79 


156% 


HDF-S- 1 


H 


B 


0.77 


0.187 


21.47 


103% 


HDF-S- 1 


Ks 


A 


0.82 


0.212 


21.34 


208% 


HDF-S-2 


.1 


A 


0.83 


0.159 


21.98 


88% 


HDF-S-2 


H 


A 


0.84 


0.181 


21.69 


112% 


HDF-S-2 


Ks 


A 


1.08 


0.176 


21.20 


198% 


HDF-S-3 


.1 


A 


1.17 


0.060 


22.07 


33% 


HDF-S-3 


K s 


B 


0.79 


0.127 


20.80 


67% 



Table 9. Main characteristics of the catalogs extracted from the final images. 



Field Passband Eff. area # Objects ra/,„, (Vega) 



CDF-S- 1 


J 


26.3 


448 


21.56 


CDF-S- 1 


Ks 


25.9 


362 


19.56 


CDF-S-2 


J 


25.6 


453 


21.78 


CDF-S-2 


Ks 


25.2 


320 


19.67 


CDF-S-3 


J 


26.3 


459 


21.45 


CDF-S-3 


Ks 


26.3 


352 


19.76 


CDF-S-4 


J 


25.9 


431 


21.68 


CDF-S-4 


H 


25.9 


227 


19.95 


CDF-S-4 


Ks 


25.9 


230 


19.13 


HDF-S- 1 


J 


25.9 


460 


21.59 


HDF-S- 1 


H 


25.9 


442 


19.64 


HDF-S- 1 


K s 


25.9 


496 


19.30 


HDF-S-2 


.1 


26.3 


672 


21.87 


HDF-S-2 


H 


25.6 


564 


19.81 


HDF-S-2 


Ks 


26.3 


521 


20.10 


HDF-S-3 


.1 


25.2 


262 


20.81 


HDF-S-3 


K s 


26.3 


395 


19.34 



scribed in Paper I, which also describes the catalog format. 
The final science-grade catalogs presented include only objects 
with S/N > 5 as computed from the magnitude error. Table [9] 
summarizes the characteristics of the extracted catalogs. The 
table gives in Col. 1 the field name; in Col. 2 the passband; in 
Col. 3 the effective area in square arcmin; in Col. 4 the num- 
ber of objects; and in Col. 5 the 80% completeness limit in the 
Vega system, mn m , as estimated from mock catalogs. 

Each catalog and its associated product log, including ex- 
tensive comparison with other data previously ingested in the 
associated database of the EIS system, takes about 2.5 min, 
corresponding to an overall mean data rate of 0.03 Mpix/s. 

5. Discussion 

The EIS-DEEP infrared data presented here was used as the 
test case for the implementation of both the stand-alone version 
of the image processing library EIS/MVM, and the EIS Data 



Reduction System as a whole. The latter is responsible for a 
variety of automatic and un-supervised processes ranging from 
the creation of the reduction blocks and of the reduced images 
on a nightly basis to the creation of stacks and their associ- 
ated final science-grade catalogs presented in the previous sec- 
tion. Because the EIS-DEEP infrared dataset is at the same time 
relatively small, enabling extensive tests to be done, and large 
enough to sample a variety of problems that emerge under real 
observing conditions, thereby helping identify trouble-spots in 
the code, it provides an excellent test case. Furthermore, this 
dataset has been reduced b oth by the EIS team, using differ- 
ent so ftware packages (Ida Costa et all Il998t iReneehnk et all 
1998) or versions of the same software JVandame et all l200lL 
and the present paper), and independently by other groups us- 
ing different techniques, thereby allowing the comparison of 
different versions. Finally, several other datasets covering the 
regions considered have been compiled since the completion 
of the EIS-DEEP observations, using either the same or differ- 
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Fig. 3. Comparison between the measured J (upper panel), and 
K s (l ower panel) magnitude s of objects in common with those 
from Vandam e et al.l J200ll) . The large filled circles represent 
the mean magnitude difference computed in bins of 1 mag, 
while the error bars denote the scatter in each bin. All mag- 
nitudes are in the Vega system. 



Fig. 2. Mosaic image of the CDF-S in the J (upper) and K s - 
band (lower). 



ent instruments with different limiting magnitudes. The avail- 
ability of these various datasets discussed below is in marked 
contrast to the time these data were first released by the EIS 
team. 

The availability of this large array of different datasets cov- 
ering the same region allows a careful assessment of the per- 
formance of the EIS system and the quality of its final science 
products by direct comparison of the position, flux scale and 
absolute calibration of objects in common between source cat- 
alogs extracted from different final images. Additionally, sta- 
tistical measures, here limited to number counts of galaxies 
and stars, compared to those computed by other authors or pre- 
dicted by models are used to assess the overall consistency with 
results from the literature. 

5. 1. Comparison with other reductions 

The comparison with other reductions of the same dataset has 
been extensively used to show, a posteriori, the shortcom- 
ings of previous releases of this dataset either produced by 
the original implemen tation of the program Jitter of Eclipse 
foevillard et all Il999l) or by an earlier version of the new 
EIS/MVM image processing library. Below, the same strat- 



egy is adopted to evaluate the performance of the current ver- 
sion of the system also used in the reduction of the infrared 
part of the DPS surve y to be presented in a fo rthcoming pa- 
per ( OlsenTTall 120061) and of the optical data in lDietrich et alJ 
(l200.4 : lMignano et alJfcOOfifc . 

As a starting point, Fig. [3] shows the difference in magni- 
tude of objects identified in images produced u sing an earlier 
versio n of the EIS/MVM library, as reported in Van dame et alJ 
d200ll) . and those of the most recent version of the package, in 
all four fields of the CDF-S region. The objects were matched 
by positional coincidence, using a search radius of 1 arcsec, 
and only 5o detections are considered. These values are used 
throughout this paper. The large filled circles are the mean 
value of the magnitude differences in bins of one magnitude, 
and the error bars indicate their scatter. The vertical dashed 
lines in the figure denote the brightest magnitude m« m , of the 
four fields considered, where mn, n is the magnitude at which 
the completeness of the catalog drops to 80%. From the top 
panel one finds that the error in the earlier version of the image 
processing library had little effect in the relative photometry 
of the /-band data, with the possible exception of the faintest 
magnitude bin, which shows the hallmark of lost flux in the 
older reduction. Surprisingly, after re-calibration one also finds 
that the J magnitudes for one of the fields exhibits an offset 
with respect to the others of as much as 0.1 mag. Closer in- 
spection of the data shows that this is associated to the field 
CDF-S-4. From the bottom panel, one finds that in the ^Ty-band 
the effects of flux loss are slightly larger but still of the order of 
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Fig. 4. Difference in magnitudes J (top panels) and K s (lower panels) as measured reducing the infrared data using the DIMSUM 
package (Daddi private communication) and the EIS software. As before, the large filled circles are the mean values of the 
difference in bins of one magnitude and the error bars are the scatter. 



< 0.05 mag, significantly smaller than the originally reported 
bias, which reached 0.15 mag at the faint end. 

To further verify the quality of the produced stacks, the J- 
and /r s -band images of the CDF-S fields 1 and 2 reduced using 
the IRAF package DIMSUM (Daddi, private communication) 
were used to extract catalogs, in the same way as those created 
by the EIS system. The catalogs were matched and the magni- 
tudes of the objects in common compared. In this comparison, 
the absolute calibration has been neglected and any residual 
offset at the bright end between the two sets of catalogs was ar- 
bitrarily removed. The results are shown in Fig.|4]with the top 
panels showing the results for the /-band and the bottom panels 
for the /T s -band, for the fields indicated in each panel. As can be 
seen from the comparison in both bands and fields the results 
are excellent, demonstrating that the earlier problems identified 
in the image processing system have been properly addressed. 

Similar results were obtained comparing other reductions 
of SOFI and ISAAC data, reduced using either the same 
or different software (e.g. Saracco, private communication; 
Dickinson, private communication). In summary, it has been 
shown that the current reductions of the infrared data obtained 
using the EIS Data Reduction System is free from any flux 
bias at faint magnitudes, yielding results in excellent agree- 
ment with those obtained using other software packages de- 
signed specifically to deal with infrared data. 



5.2. Comparison with other datasets 

The validation of the infrared reductions presented here 
can be extended, comparing them to other fully calibrated 
datasets available in the literature and covering the same re- 
gions surveyed by EIS-DEER Examples of these datasets in- 
clude: 1) data from the Two Mass All Sky Survey (2MASS, 
iKleinmann et all 1 19941) in all three passbands, useful to check 
the bright end; 2) de ep data taken with ISAAC at the VLT 
for the FIRES survey ILabbe et all 12003 ) covering one of the 
pointings of the HD F-S region ; 3) da ta taken in H-band of 
the CDF-S region bv lMov et all J2003I): and data from t he Las 
Campanas Infrared Survey (LCIRS, IChen et all 120021) of the 
HDF-S region. 

The results of these various comparisons are shown in 
Figs.|5l|8] Fig.|5]shows the offset in position (left panels) and 
magnitude (right panels) for objects in common with 2MASS 
within 1 arcsec of EIS-DEEP detections. Each row shows the 
results obtained combining all the fields for the different bands 
with J on top and K s on the bottom. While the number of ob- 
jects is relatively small, considering the magnitude interval of 
the two samples, the main advantage of 2MASS is that data are 
available for all fields allowing the use of a single reference to 
check both the astrometric and photometric calibration. From 
Fig-HJone finds that there is a relative offset between the refer- 
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ence system of the two datasets. However, this offset is nearly 
the same for all passbands ranging from 0.27 to 0.43 arcsec 
in right ascension and from 0.19 to 0.26 arcsec in declination 
with a scatter of < 0.40 arcsec in all bands. Since the objects in 
common are predominantly from HDF-S, the observed offset 

reflects the difference between the Tycho and GSC2 coordi- 

J 1 1 ~i 

nates in the HDF-S region. As shown by (Olsen et al., 2006), 
no systematic shifts are found relative to the USNO-B catalog 
used to calibrate the different regions covered by the DPS in- 
frared dataset. 

Comparison of the magnitudes yields offsets of -0.02, 
0.004, and -0.05 mag in J, H and K s , respectively, if one dis- 
cards very bright objects which tend to be partly saturated 
in the deep data. Another important fact is that, except for 
the Malmquist bias seen in the plots, there is no evidence of 
systematic trends with magnitude in the interval considered. 
Furthermore, the small scatter suggests some degree of consis- 
tency between the calibration of the images involving as much 
as 7 fields distributed over two distant regions of the sky, in 



some cases observed in different nights. All this evidence sug- 
gests a consistent absolute calibration of the EIS-DEEP data. 

The behavior at faint magnitudes is evaluated by comparing 
the present data to the deep observations of the FIRES survey 
(Labbe et al., 120031) . Since the authors do not provide the data 
with an astrometric solution a crude shift has been applied to 
obtain a proper match between the final images derived for the 
two datasets. Therefore, in this case only magnitudes are com- 
pared. Furthermore, in the comparison only objects detected at 
5o in both catalogs are considered. Fig.|6]shows the result of 
this comparison for all the available passbands. From the figure 
one finds that regardless of the passband, there is no systematic 
deviation of the magnitude differences as a function of the mag- 
nitude apart from a constant offset. One possible exception is 
the difference at the faintest bin in the /-band, which seems to 
reflect the effects of a Malmquist bias rather than flux-loss at 
the faint end. The comparison yields offsets of -0.1 1, 0.06 and 
0.001 mag in /, H and K s , respectively. 
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Fig. 8. Same as Fig.[7]but now comparing positions and magnitudes measured by the LCIRS of the HDF-S. 



In the //-band t wo more datasets are available for compar- 
ison. One is that of Mo v et alJ (|2P03) covering the CDF-S re- 
gion, for which there is an overlap with the somewhat shallow 
EIS-DEEP observations of CDF-S-4, and the other is that ob- 
tained as part of the LCIRS co verage of th e HDF-S region. The 
results of the comparison with Mov et alJ d2003l) are shown in 
Fig.0 where the differences in position (left panel) and magni- 
tude (right panel) of objects in common are plotted. The mean 
offset in position is less than 0. 1 arcsec with an rms of about 
0.45 arcsec, consistent with both having an error of the or- 
der of ~ 0.25 arcsec. This result shows that the positions in 
the CDF-S region have no systematic shift. The magnitudes 
are also in good agreement except for an offset of 0.1 mag, 
somewhat larger than the 0.06 mag derived from the compar- 
ison with FIRES data of HDF-S. The origin of this offset is 
unknown and may be due to an inadequate absolute calibra- 
tion. However, this calibration can only be improved by addi- 
tional observations under photometric conditions and a better 
calibration plan. Again there is no systematic deviation of the 
magnitude differences over the entire magnitude interval con- 
sidered. This is in marked contrast with the results obtained 
comparing the EIS data with the LCIRS //-band data taken in 
the region of the HDF-S, shown in Fig. [8] Here the magnitude 
difference tends to increase for fainter objects, with the LCIRS 
magnitudes being fainter than those assigned by the EIS sys- 
tem. This trend reaches about 0.2 mag relative to the bright end, 
except for the faintest bin shown affected by the Malmquist 
bias. Since this effect is not seen in any of the other compar- 
isons discussed above, it suggests that this effect stems from 
the LCIRS reductions. Note that the offset in the bright end is 
0.12 ± 0. 07, consistent with that computed from the compar- 
ison with Mo v et alJ (120031) of 0.1 mag, but somewhat larger 
than the offset of 0.06 mag measured against FIRES. Still, re- 
gardless of the exact value this result also suggests a remark- 
able consistency in the photometric calibration of the CDF-S 



and HDF-S regions. On the other hand, the left panel of Fig. [5] 
shows the same offset in the coordinates as observed relative 
to the 2MASS data. Note that the absolute differences can eas- 
ily be corrected. To construct optical/infrared catalogs the as- 
trometric calibration of infrared images should be done using 
catalogs extracted from the optical images, thereby minimizing 
possible relative offsets. 

5.3. Comparison of galaxy and star counts 

It is not always possible to assess the quality of the data by di- 
rect comparison with overlapping observations. An alternative 
is to compute different statistical measures, the simplest exam- 
ples being galaxy and star counts, and compare them with those 
obtained by other authors in other regions of the sky, typically 
used for galaxies, or for stars with those predicted by models 
of the Galaxy. 

Figs.|9]and 1 101 show the galaxy counts computed from the 
catalogs extracted from the EIS-DEEP images of each field 
(columns) and passband (rows) of the HDF-S and CDF-S re- 
gions, respectively. This is done in order to minimize differ- 
ences that could be introduced in mean counts from field- 
to-field variations in depth of the observations. Here all ob- 
jects with stellarity index < 0.95 were classified as galax- 
ies. The present counts are compared to those published ear- 
lier by different authors ( Vaisanen et al., 2000; Martini, 12001 [ 
llovino etail l2005h . The good agreement between the various 
counts is further evidence of the scientific quality of the data 
and calibration presented in this paper. One can also see the 
strong field-to-field variations in limiting magnitude, a problem 
one must take into account when dealing with mosaics cover- 
ing large areas of the sky. 

For completeness, Figs. ^2 an d Elshow similar plots for 
the stellar counts, comparin g them wi th the predictions of the 
Galactic model o nGirardi et al. (2005). Considering the uncer- 
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Fig. 9. Galaxy number counts computed from the present data 
for each band (rows) and field (columns) of the HDF-S re- 
gion (solid circles) compar ed to the publi s hed v alues with 
diamonds denoting those of Vaisanen et al. (2000); triangles 
iMartinil J200ll) and sauares llovino et al.ld2005l) ~ 



Fig. 11. Star counts in the three passbands (rows) for the 
three fields (columns) covering HDF-S (full circles) compared 
to the predictions (solid line) computed from the model of 
lGirardietalJll2005l) . 



tainties involved in the model and the shot-noise the agree- 
ment is quite remarkable, with the possible exception of the 
K s counts in the CDF-S-1 and -3. 
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Fig. 10. Same as Fig.|9]but for the CDF-S region. 
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Fig. 12. Same as Fig.[TT]but for the four fields of CDF-S 

6. Summary 

This paper presents the results of the infrared part of the EIS- 
DEEP survey carried out in the period 1998-1999 covering the 
HDF-S and CDF-S regions. The paper includes re-processed 
and re-calibrated images of previously released data, as well 
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as new sets of data not released before. These data are avail- 
able at the CD S 2 and supersede earl i er releases descr i bed in 
the preprints by da Costa et al. ( 199^: lRengelink et alJ fl99^) 
and lVandame et alJ ( 1200 ll) which represent different phases of 
the project. In contrast to the original release, new datasets have 
become available, allowing a critical examination of the results 
by carrying out an array of direct comparisons with reductions 
of the same dataset as well as comparison with independent 
datasets. 

The results show that the past problems with flux loss of 
Jitter and an earlier version of the EIS/MVM library have been 
successfully corrected with the magnitudes measured from the 
EIS images showing no systematic deviations relative to those 
derived from images reduced using other reduction techniques. 
Furthermore, comparisons with data obtained by other authors 
also suggest that the absolute calibration is reliable to within 
~ 5 — 10%. To obtain a better accuracy would require addi- 
tional observations under photometric conditions and a good 
sampling of standards over the entire night. 

The results obtained in this paper, combined with those dis- 
cussed in previous papers of this series, demonstrate the power 
and versatility of the EIS Data Reduction System which can 
be used to efficiently process optical and infrared data in the 
same environment yielding reliable results. The EIS-DEEP sur- 
vey was the precursor of the the Deep Public Survey (DPS), 
which combines optical wide-field imaging with infrared data. 
The results of this survey will be reported in forthcoming pa - 
pers of this series Joisen et all l2006t iMignano et all EoOfil) . 
Combined, these data provide information in up to seven bands 
from the ultraviolet to the near-infrared, thus being an excel- 
lent test case for future optical/infrared surveys. The results 
presented here and in other papers of this series demonstrate 
that the EIS Data Reduction System meets the requirements 
for long-term programs using the VST and VISTA survey tele- 
scopes. 
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